Interface magnetism in Fe203/FeTi03-heterostructures 
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To resolve the microscopic origin of magnetism in the Fe203/FeTi03-system, we have performed 
density functional theory calculations taking into account on-site Coulomb repulsion. By varying 
systematically the concentration, distribution and charge state of Ti in a hematite host, we compile 
a phase diagram of the stability with respect to the end members and find a clear preference to 
form layered arrangements as opposed to solid solutions. The charge mismatch at the interface 
is accommodated through Ti*^ and a disproportionation in the Fe contact layer into Fe^"'', Fe^^, 
leading to uncompensated moments in the contact layer and giving first theoretical evidence for the 
lamellar magnetism hypothesis. This interface magnetism is associated with impurity levels in the 
band gap showing halfmetallic behavior and making Fe2 03/FeTi03 heterostructures prospective 
materials for spintronics applications. 



A challenge of today's materials science is to 
design ferromagnetic semiconductors operating at 
room-temperature (RT) for spintronics devices. 
Most of the efforts concentrate on homogeneous 
doping of semiconductors with magnetic impuri- 
ties [l[ H, H, II] , but the interfaces in complex oxides 
prove to be another source of novel behavior [1,[1,[3|- 
The unique magnetic properties of the hematite- 
ilmenite system (a canted antiferromag- 

net and a RT paramagnet, respectively) currently 
receive revived interest as a possible cause of mag- 
netic anomalies in the Earth's deep crust and on 
other planets jll| as well as for future device appli- 
cations [iElIai. 

Both hematite (a = 5.035A, c = 13.75lA (H) 
and ilmcnite (a = 5.177A, c = 14.265A [ll]) crys- 
tallize in a corundum(-derived) structure shown in 
Fig.[T]where the oxygen ions form a distorted hexag- 
onal close packed lattice and the cations occupy 
2/3 of the octahedral sites. In a-Fe203 (space 
group i?3c) there is a natural modulation of elec- 
tronic density along the [0001]-direction where neg- 
atively charged 30^^ layers alternate with positively 
charged 2Fe"^+ layers. At RT the magnetic moments 
of subsequent iron layers couple antiferromagneti- 
cally (AFM) in-plane with a small s pin- canting, at- 
tributed to spin-orbit couphng Il7l. [isj . In il- 
menite, FeTiOs, Fe- and Ti-layers alternate, reduc- 
ing the symmetry to i?3 and the corresponding se- 
quence is 302-/2Fe2+/302-/2Ti'*+ with AFM cou- 
pling between the Fe layers and Tn =56-59 K (Toj . 

At the interfaces (IFs) in hematite-ilmenite ex- 
solutions, charge neutrality is disrupted. One way 
to balance the excess charge at the interface is by 
a disproportionation in the Fe layer, now becom- 



ing mixed Fe^+ and Fe'^"'". This lamellar mag- 
netism hypothesis (LMH) was proposed by Robin- 
son et al. [lOl based on bond valence models and 
kinetic Monte Carlo simulations (kMC) with empir- 
ical chemical and magnetic interaction parameters. 
The increased technological interest in this system 
calls for an atomistic material specific understanding 
that can only be obtained from first principles calcu- 
lations. A previous density functional theory (DFT) 
study within the generalized gradient approximation 
(GGA) found no evidence for the LMH [2l|. How- 
ever, electronic correlations, not included in the lo- 
cal (spin-) density approximation (LSDA) or GGA 
of DFT, play an important role in transition metal 
oxides. Such effects were considered recently within 
LSDA-hU [H or using hybrid functionals [Hi for 
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FIG. 1: (Color online) Crystal structure of FeTi03, 
showing half of the 60 atom unit cell. The cation sites are 
numbered, oxygen occupies the edges of the octahedra. 



single Ti impurities in hematite, however layered ar- 
rangements and interfaces were not addressed. 

In this paper we have performed DFT calculations 
including a Hubbard U [13] for the end members 
Fe203 and FeTiOs, as well as their interfaces and 
solid solutions (SS). By varying systematically the 
concentration, distribution and charge state of Ti 
incorporated in a Q;-Fe2 03-host, we explore different 
scenarios for the charge compensation mechanism 
and its consequences for the magnetic and electronic 
behavior. Finally, we compile a phase diagram of the 
stability of the different configurations with respect 
to the end members as a function of Ti-doping also 
taking into account the effect of strain. 

Our DFT-GGA (l^ calculations are performed 
using the all-electron full-potential augmented plane 
waves (FP-APW) method as implemented in 
WIEN2k [25]. Electronic correlations are considered 
within the fully locahzed hmit (LDA-hU) [24]. The 
systems are modeled in the hexagonal primitive unit 
cell (shown in Fig. [T]) containing 30 and 60 atoms. 
For these 24 and 15 fc-points in the irreducible part of 
the Brillouin zone were used, respectively. Inside the 
muffin tins (i?^^"^' = 1.80 bohr, i?^T ~ 1.60 bohr) 
wave functions are expanded in spherical harmonics 
up to Z^ax — 10 and non-spherical contributions to 
the electron density and potential up to lf°^^ = 6 
are used. The energy cutoff for the plane wave rep- 
resentation in the interstitial is -E^ax = 19 for 
the wave functions and i?J^ax = 196 Ry for the po- 
tential. These convergence parameters ensure a nu- 
merical accuracy of energy differences better than 
O.OleV/60 atom cell. A full structural optimization 
of internal parameters has been performed (27j . 

As a starting point we have modeled the end mem- 
bers Fe203 and FeTiOa. In agreement with previous 
calculations [13, HI], GGA-I-U considerably improves 
the band gap of hematite from 0.43 eV (GGA) to 
2.2 eV for U=6 eV and J=l eV in close agreement 
with measured values of 2.14-2.36 eV [2^, [13]. Also 
the type of band gap changes from a Mott-Hubbard 
between Fe3c?- Fe3d states to charge transfer after 
the scheme of Zaanen et al. [3l| between occupied 
02p and empty FeSd-states. 

For ilmenite GGA incorrectly predicts a metal- 
lic state jsl], hence the inclusion of Hubbard U 
{U = 8 eV, J = 1 eV) is decisive to obtain a 
Mott-Hubbard gap of 2.18 eV (A'^^p = 2.54 eV [9]) 
between the occupied Fe(i22 -orbital and the unoc- 
cupied Fe3d states in one spin channel (all Fe3d 
orbitals being occupied in the other spin-channel). 
These U and J values are used both on Fe and Ti in 
the following. A Fe^+ZTi'^+charge arrangement lies 
0.63 eV/p.f.u. above the ground state Fe'^+/Ti'^+. 
For Fe^+Ti'^+Oa the AFM and FM couphng between 
Fe-layers is nearly degenerate, consistent with the 



TABLE I: Relative stability (eV/60 at. cell) of the differ- 
ent cation arrangements in the Ti**^, Fe'^^, Fe^"'' charge 
state (strained at the Fe2 03 lattice parameters) with re- 
spect to the most stable configuration whose energy is 
set to 0.0 eV. The positions occupied by Ti are denoted 
as subscripts according to Fig. [1] (SSL: spin-sublattice) . 
The total magnetic moment, electronic behavior (hm/m 
denotes halfmetallic/metallic] as well as Ti*"*"-© and 
Fe^"*'-0 distances at the interface are also displayed. 
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low Tat. 

In the following we vary the concentration and 
distribution of Ti in a Fe203-host. The positions 
of the Ti-ions are given as subscripts and follow the 
notation of Fig. (TJ e.g. T3_4^7^g^ii_i2 describes pure 
ilmenite. Table [J contains the energetic stability, 
structural, magnetic and electronic properties of dif- 
ferent cation arrangements and concentrations. We 
start the discussion with Ilmaa which corresponds to 
four Ti-ions out of 24 cations in the 60-atom cell. 
We find that the formation of a compact ilmenite- 
like block with a Fe-layer sandwiched between two 
Ti-layers (T3^4), is by 0.36 eV more favorable than 
incorporation of single Ti layers in the hematite host 
(Ti,2,i3,i4)- The spin density of T3_4 plotted in 
Fig. [2^) shows that the central Fe-layer turns into 
Fe^+ and the charge mismatch at the IF is compen- 
sated by Fe^+jFe'^"'" in the contact layer, giving theo- 
retical evidence from first principles for the lamellar 
magnetism hypothesis of Robinson et al. [2^. Our 
GGA-fU calculations show that Ti*+4 shares faces 
with Fe^+5, while Fe^+g shares faces with Fe'^+s 
from the next hematite layer. Such a configura- 
tion was proposed using bond-valence sums [33| and 
kMC [33] only after considering both chemical and 
magnetic interactions. 

The formation of layered arrangements (T3.4) is 
favored compared to a more random distributions 
with 50% substituted cation layers (e.g. T^^iq, Ts^y 
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FIG. 2: (Color online) a) Spin density distribution and 
b-e) total and projected density of states of Ti double 
layer in hematite (T3,4) with Fe'^'''i,2, Ti''+3,4 and a dis- 
proportionated Fe^+s, Fe'^+e layer at the interface. The 
position of the Fermi level is set to 0.0 eV and denoted by 
a short line. In a) positions of Ti, Fe and O are marked 
by light grey, red and small black circles, respectively. 



or Ts^s)- With respect to magnetism, each Ti 
ion adds a magnetic moment of independent 
of whether the extra electron is localized at Ti 
(Ti^+) or at a neighboring Fe (Fe^+). In solid so- 
lutions the total magnetic moment depends on the 
site and sublattice where Ti is built in. We find 
that incorporation in the same spin-sublattice Ts^io 
(which maximizes the magnetic moment) is favored 
by 0.22 eV compared to the AFM T5 j. Taking into 
account local lattice relaxations enhances the energy 
gain compared to previous calculations by Velev et 
al. (0.08 eV)[23]. Still, some degree of Ti disorder 
is likely in quickly cooled samples, reducing the ex- 
pected magnetization as observed by Chambers et 
al. 3 (0.5±0.15/Xij/Tifora;Ti = 0.15). For higher Ti 
concentrations and longer annealing steps the forma- 
tion of the thermodynamically more stable layered 
ferrimagnetic phase is expected, consistent with the 
strong correlation between cation order and ferri- 
magnetism found in annealed samples [1, and 
the saturation magnetization of S/is/mol measured 
in epitaxial films with xti < 0.63 [12]. While be- 
low the ordering temperature of ilmenite only il- 
menite lamella with an odd number of Ti-layers are 
expected to carry a non-zero magnetization, above 
56 K the net magnetic moment will be solely due to 
the uncompensated magnetic moments in the con- 
tact layers, independent of the number of Ti-layers 
within the paramagnetic ilmenite lamella. 

Concerning the electronic properties, doping 
Fe203 with Ti leads to impurity levels in the 
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FIG. 3; (Color online) Formation energy (eV/p.f.u) with 
respect to the end members as a function of ilmenite 
concentration x. Layered arrangements (solid solutions) 
are denoted by full (open) symbols. 



band gap arising from the occupied (i^2 -orbital of 
Fe^+ ions in the contact layer. The density of 
states plotted in Fig. [2}D-e) shows that these states 
are pinned at the Fermi level, leading to fully 
spin-polarized carriers and half-metallic behavior for 
T3^4. This trend is robust with respect to U [s^ 
and is observed for most of the studied cation con- 
centrations and arrangements after structural relax- 
ation. Experimentally, semiconducting behavior and 
a drop in resistivity of several orders w.r.t. the end 
members is measured [1, [l^, [l^l with values sug- 
gesting localized rather than itinerant carriers con- 
sistent with the picture we obtain from LDA+U. 

Next we turn to the energetic stability as a func- 
tion of the Ti-concentration displayed in Table U and 
the phase diagram in Fig. [31 Charge compensation 
through Ti^+ and disproportionation of iron into 
Fe^+ , Fe'^+ is strongly favored compared to compen- 
sation involving Ti'^"'", especially after optimization 
of the internal structural parameters. Moreover, the 
formation of layered configurations (full symbols) is 
preferred over disordered arrangements (empty sym- 
bols) except for very high (> 83%) and very low 
concentrations (< 17%), consistent the miscibility 
gap from thermodynamic data (e.g. [21I, H^). The 
linear increase of formation energy in the range be- 
tween Ilmi7 and Ilmgg indicates that straining Ti 
doped Fe203 to the Fe203 lattice parameters gets 
increasingly unfavorable with growing x. On the 
other hand, using the ilmenite lattice parameters at 
X = 66% instead of hematite (volume increase of 
8.7%) leads to an energy gain of 0.22 eV/p.f.u. for 
T7,8,ii,i2 (rsd filled up-triangles in Fig. [3]). 

An interesting trend is observed in the shortest 
cation-oxygen bond lengths (cf. Table [l|, which 
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tend to relax towards the values in the respective 
end member. While dpo3+_o (not shown) remains 
close to the value in bulk hematite (1.96 A), the bond 
lengths of the Ti-impurity and the neighboring Fe^+ 
relax towards the values in bulk ilmenite (1.92 and 
2.O7A, respectively). 

In Fe doped ilmenite the trend towards layered 
arrangements is retained for 66%, e.g. T7_8,ii,i2 is 
favored by 0.21 eV compared to T3_8,ii,i2i but or- 
dered and disordered phases are nearly degenerate 
at 83%. Fe substituting for Ti in the ilmenite lat- 
tice is Fe'^"'' . Additionally, one Fe in the neighboring 
layer turns Fe'^^ to compensate the charge, forming 
a Fe^+, Fe'^+contact layer. The substituted Fe shows 
a very strong tendency to couple antiparallcl to the 
neighboring Fe-layers. 

In summary, we present a comprehensive 
GGA-t-U-study of the cation, charge and magnetic 
order in the hematite-ilmenite system showing a 
strong preference towards formation of layered con- 
figurations as opposed to solid solutions. At the 
interface between hematite and ilmenite blocks we 
find evidence for the lamellar magnetism hypothe- 
sis [13] with a disproportionated Fe^"*", Fe'^+contact 



layer to accommodate the polar discontinuity. These 
uncompensated moments lead to ferrimagnetic be- 
havior of the system. The (1^2 orbital in one spin- 
channel (all d orbitals being occupied in the other) at 
the Fe^+-sites in the contact layer crosses the Fermi 
level leading to halfmetallic behavior in most of the 
studied compositions. The ferrimagnetism emerg- 
ing at the interface of two antifcrromagnetic oxides 
like hematite and ilmenite is an impressive example 
of the novel functionality that can arise as a conse- 
quence of a polar discontinuity and the richer pos- 
sibilities to compensate it that complex oxides of- 
fer. Recently, an exchange bias of more than 1 Tesla 
was reported in this system [37]. Further phenom- 
ena such as oscillatory exchange coupling and spin- 
polarized transport remain to be explored in con- 
trolled epitaxial Fe203-FeTi03-multilayers on the 
route to possible device applications. 
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